Magnetization of a strongly interacting two-dimensional electron system in 

perpendicular magnetic fields 
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We measure the thermodynamic magnetization of a low-disordered, strongly correlated two- 
dimensional electron system in silicon in perpendicular magnetic fields. A new, parameter-free 
method is used to directly determine the spectrum characteristics (Lande g-factor and the cyclotron 
mass) when the Fermi level lies outside the spectral gaps and the inter-level interactions between 
quasiparticles are avoided. Intra-level interactions are found to strongly modify the magnetization, 
without affecting the determined g* and m* . 
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Magnetization is one of the least studied properties 
of two-dimensional (2D) electron systems: signals asso- 
ciated with the magnetization of 2D electrons are weak, 
and measuring them is a challenging experiment. Few ex- 
perimental observations of the de Haas- van Alphen effect 
in 2D electron systems were made using SQUID magne- 
tometers pj, pick up coils lithographed above the gate 
Q, or torque magnetometers A novel method has 
recently been used by Prus et al. Q and Shashkin et al. 
[f| to measure the spin magnetization of 2D electrons 
in silicon metal-oxide-semiconductor field-effect transis- 
tors (MOSFETs). This method entails modulating the 
magnetic field with an auxiliary coil and measuring the 
imaginary (out-of-phase) component of the AC current 
induced between the gate and the 2D electron system, 
which is proportional to d\ijdB (where /z is the chem- 
ical potential). Using the Maxwell relation, d[ijdB = 
—dM/dn s , one can then obtain the magnetization M by 
integrating the induced current over the electron density, 
n s . Pauli spin susceptibility has been observed to be- 
have critically near the 2D metal-insulator transition, in 
agreement with previous transport measurements ja, |j| . 

Here we apply a similar method to study the ther- 
modynamic magnetization of a low-disordered, strongly 
correlated 2D electron system in silicon MOSFETs in 
perpendicular and tilted magnetic fields. By measur- 
ing dfi/dB at non-integer filling factors, we directly de- 
termine the spectrum characteristics without any fitting 
procedures or parameters. As compared to previously 
used measuring techniques, the remarkable advantage of 
the novel method is that it probes the spectrum of the 
2D electron system with the Fermi level lying outside the 
spectral gaps so that the effects of interactions between 
quasiparticles belonging to different energy levels (inter- 
level interactions) are avoided. Although intra-level in- 
teractions are found to strongly affect the magnetization, 
the extracted Lande g-factor and the cyclotron mass are 
insensitive to them. Therefore, measured spectrum char- 
acteristics are likely to be identical with those of a contin- 
uous spectrum. The so-obtained ^-factor has been found 



to be weakly enhanced and practically independent of 
the electron density down to the lowest densities reached 
(« 1.5 x 10 11 cm -2 ), while the cyclotron mass becomes 
strongly enhanced at low n s . 

Measurements were made in an Oxford dilution refrig- 
erator on clean (lOO)-silicon samples with peak electron 
mobilities of 3 m 2 /Vs at 0.1 K and oxide thickness of 
149 nm. Magnetic field B was modulated with a small 
AC field B mod in the range of 0.005 - 0.03 T at a low 
frequency / = 0.05 — 0.45 Hz to minimize possible me- 
chanical resonances and avoid overheating the sample. 
The latter was verified by monitoring the temperature- 
dependent sample resistance at B = and in the range 
of filling factors 1 to 6 with and without modulation. 
Noticeably higher amplitudes B mo( \ and/or frequencies 
/ caused overheating of the mixing chamber and were 
avoided. The in-phase and out-of-phase components of 
the current between the gate and the 2D electron system 
were measured with high precision (~ 10~ 16 A) using a 
current- voltage converter and a lock- in amplifier. The 
imaginary (out-of-phase) current component is equal to 
Im i = (2TrfCB mo< i/e)dfj,/dB, where C is the capaci- 
tance of the sample. For measurements of the capaci- 
tance, a similar circuit was used with a distinction that 
the gate voltage was modulated and thus the imaginary 
current component was proportional to the capacitance. 
The electron density was determined from the capaci- 
tance oscillations. 

Typical experimental traces of the gate current in a 
perpendicular magnetic field of 5 T are displayed in 
Fig. ^ Sharp dips in the out-of-phase component, seen 
at integer filling factors v = n s hc/eB±, reflect gaps in 
the density of states: dips at odd filling factors corre- 
spond to the valley splitting, the ones at v = 2 and 6 
are due to the spin splitting, and the dip at v — 4 is due 
to the cyclotron splitting. However, there are no cor- 
responding features in the in-phase current component, 
which ensures that we reach the low-frequency limit and 
the measured d[ijdB is not distorted by lateral trans- 
port effects. This is further confirmed by the fact that 
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FIG. 1: Out-of-phase (solid line) and in-phase (dotted line) 
current components as a function of the electron density 
in a perpendicular magnetic field of 5 T and T — 0.8 K. 
-Bmod = 0.022 T and / = 0.45 Hz. The value dfi/dB is indi- 
cated in units of the Bohr magneton fiB- In the right-hand 
inset, we demonstrate proportionality of Im i to frequency: 
the solid and dashed lines (vertically shifted for clarity) cor- 
respond to 0.45 and 0.1 Hz, respectively; the {/-component of 
the latter is multiplied by 4.5. The left-hand inset illustrates 
magnetization per electron. 

the out-of-phase current is proportional to the excitation 
frequency as displayed in the right-hand inset to Fig. ^ 
Magnetization per electron can be extracted by integrat- 
ing the measured out-of-phase signal with respect to n s , 
as shown in the left-hand inset to Fig. ^ for illustration. 
The magnetization exhibits the expected sawtooth oscil- 
lations, with sharp jumps at integer filling factors (note 
that the height of the jumps yields values that are smaller 
than the level splitting by the level width). 

If the disorder and interactions are disregarded, in 
quantizing magnetic fields (except at integer filling fac- 
tors) the derivative dfi/dB — —8M/dn s is equal to 

dfi |Y 1 . \ 2m e ,11 

where /is is the Bohr magneton, N is the Landau level 
number, m e and mi, = 0.19m e are the free electron mass 
and band mass, respectively, and go = 2 is the (/-factor in 
bulk silicon. Disorder smears out the dependences which 
otherwise would consist of a series of delta- functions. In- 
teractions modify this picture in two ways: (i) by renor- 
malizing the values of the cyclotron mass and (/-factor 
and (ii) by providing a negative contribution of order 
— {e 2 /eIb)^} 1 / 2 to the chemical potential H, || (here 
e is the dielectric constant, Is is the magnetic length, 
and {v} is the deviation of the filling factor from the 
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FIG. 2: (a) Capacitance in B — 8 tesla and in B — as 
indicated. The (noise averaged) geometric capacitance is de- 
picted by a dashed line, (b) Im i oc d/i/dB in a perpendicular 
magnetic field of 8 tesla. The maximum values possible in 
a non-interacting system (see text) are depicted by dashed 
lines. 

nearest integer). The latter effect, which is caused by 
the intra-level interactions between quasiparticlcs, leads 
to the so-called negative thermodynamic compressibility 
near integer filling factors predicted by Efros |8( and ex- 
perimentally observed in Refs. [lfllllj . 

In Fig. |2 we compare capacitance C with d/i/dB, 
measured at the same magnetic field value and plotted 
versus n s around the filling factor v = 2. The capac- 
itance consists of two contributions: 1/C = l/C gco + 
l_lAe 2 (dn s /dfi), where C goo is the geometric capacitance 
|12| depicted by the dashed line in Fig. |21 (a), and A is 
the sample area. (Note that the geometric capacitance 
slightly increases with n s since the thickness of the 2D 
electron layer — and, therefore, the average distance be- 
tween the 2D layer and the gate — decreases with the 
gate voltage.) The second term is responsible for the 
dip centered at n s — 3.87 x 10 11 cm~ 2 , corresponding to 
v = 2, and sharp maxima on both sides of it. Note that at 
these maxima, the capacitance exceeds C' gco , which cor- 
responds to the negative thermodynamic compressibil- 
ity discussed above. Farther from integer filling factors, 
the intra-level interaction corrections become weak, being 
proportional to {j/}^ 1 / 2 , and the measured capacitance 
approaches C gco (as long as the broadening of Landau 
levels is negligible, i.e., dn s /d^, ^> dn s /dn\ B=Q ). 

Similar maxima on both sides of v — 2 are seen in 
the magnetization data shown in Fig. [21 (b). At the 
maxima, the derivative d/i/dB exceeds maximum val- 
ues possible in a non-interacting 2D electron gas, which 
are determined by Eq. (1) and are depicted in the fig- 
ure by dashed lines. The possibility that d/i/dB might 
exceed its maximum non-interacting values due to intra- 
level Coulomb interactions between quasiparticles was 
predicted by MacDonald et al. in fact, this is how 
negative compressibility [j| manifests itself in magneti- 
zation measurements. Sharp spike just above v = 4 and 
maxima on both sides of v = 2 in the dependence shown 
in Fig. n are of the same nature. 

It is straightforward to obtain the effective (/-factor 
from the data for dji/dB. In accordance with Eq. (1), 
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FIG. 3: Illustration of how the effective g-factor (a, b) and the cyclotron mass (c, d) have been measured. The imaginary 
current component is plotted as a function of the deviation of the filling factor from v — 2. In perpendicular magnetic fields, 
the difference between dfi/dB for spin-down (J.) and spin-up (T) electrons yields g* in units of the Bohr magneton. In tilted 
magnetic fields, the difference between dfi/dB for electrons with N — 1 and N = is equal to 2ps (m e /m*) cos <f>. The dashed 
lines show noise-averaged values. B mo d = 0.022 T (a, b) and 0.0055 T (c, d). 



it is equal (in units of the Bohr magneton) to the differ- 
ence between d/i/dB for spin-down (J,) and spin-up (f) 
electrons belonging to the same Landau level: [Ib 9* = 
(dfi/dB)i — (dfi/dB)^. It is important that this method 
of determining the g-factor does not require the use of 
any fitting procedures or parameters. Figure [3] (a, b) 
shows measured dfi/dB as a function of the deviation 
of the filling factor from 2 at two values of magnetic 
field. Near v = 2, there are sharp intra-level interaction- 
induced structures discussed above; these regions have 
been excluded from the analysis. However, farther from 
v = 2, the dependences for v < 2 and v > 2 become par- 
allel to each other. This ensures that the so-determined 
g* is not affected by the valley splitting 0, Q and 
intra-level interaction effects H,|| discussed above. The 
latter contribute equally to both spin-up and spin-down 
dependences and cancel each other out. Disorder also 
contributes equally to dfi/dB on both sides of v = 2: we 
have found that at magnetic fields down to approximately 
3 T, there are wide regions of filling factors where capac- 
itance (i.e., the density of states) is symmetric around 
v = 2 (see, e.g., Fig. [21 (a)); furthermore, closeness of 
the capacitance to C gGO attests that the disorder-induced 
corrections are small. At lower magnetic fields, however, 
the electron-hole symmetry around v = 2 breaks down, 
which sets the lower boundary for the range of magnetic 
fields (and, consequently, electron densities). Note that 
temperature smears out the dependences in a way simi- 
lar to disorder: at higher temperatures, the capacitance 
at half- integer filling factors decreases, which leads to a 
worsening of the method accuracy. 

In Fig. 2| we plot the measured <?-factor along with the 
one previously obtained from transport measurements 
(solid line). One can see that there is no systematic de- 
pendence of the g-factor on n s : it remains approximately 
constant and close to its value in bulk silicon even at the 
lowest electron densities, which is in good agreement with 
the transport and magnetocapacitance 01 results. 

The same method can be used for determination of the 
cyclotron mass in tilted magnetic fields strong enough to 
completely polarize the electron spins [16j. If (and only 
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FIG. 4: The effective g-factor (circles) and the cyclotron 
mass (squares) as a function of the electron density. The solid 
and long-dashed lines represent, respectively, the (/-factor and 
effective mass, previously obtained from transport measure- 
ments |Te|. and the dotted line is the Pauli spin susceptibility 
obtained by magnetization measurements in parallel magnetic 
fields The critical density n c for the metal-insulator tran- 
sition is indicated. 



if) the spin splitting exceeds the cyclotron splitting, the 
gap at v = 2 lies between Landau levels 0| and If, and 
the difference (dfj,/dB)jM = i — (dfi/dB)N=o is equal to 
2/j.b (m e /m*) cos 4>, where <j) is the tilt angle. Once the 
electron spins are fully polarized at filling factors above 
v = 2, the tilt angle is automatically large enough for the 
level crossing to have occurred. The region of explorable 
electron densities is restricted from above by the condi- 
tion that the electrons must be fully spin-polarized, while 
with our current set-up, the maximum magnetic field at 
which we can apply the modulation is only 8 tesla capable 
of pola rizing the electron spins up to n* w 2 x 10 11 cm -2 
[H Il7j. Figure El (c, d) shows dfx/dB as a function of 
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\v — 2 1 under the condition n s < n* at two tilt angles 
|18|. The extracted cyclotron mass at electron densities 
1.55 and 1.35 x 10 11 cm -2 is significantly enhanced. At 
densities below 1.35 x 10 11 cm~ 2 , the symmetry of capaci- 
tance on both sides of the v = 2 gap breaks down, making 
the determination of m* impossible. As a result, we were 
only able to obtain two data points. Nevertheless, good 
agreement with the effective mass previously obtained by 
transport measurements (Fig.^J) demonstrates the appli- 
cability of the new method and adds credibility to both 
transport and magnetization results. 

We stress once again that the advantage of the new 
method we use here is that it allows determination of the 
spectrum of the 2D electron system under the condition 
that the Fermi level lies outside the spectral gaps, and 
the inter-level interactions are avoided. Being symmetric 
about v = 2, the intra-level interactions are canceled out 
in the data analysis and do not influence the extracted 
g-factor and cyclotron mass. Therefore, the obtained val- 
ues g* and m* are likely to be identical with those for 
a continuous spectrum, and the comparison with previ- 
ously found values of the g-factor and the effective mass 
is valid. 

To summarize, thermodynamic magnetization mea- 
surements in perpendicular and tilted magnetic fields al- 
low determination of the spectrum characteristics of 2D 



electron systems and show that enhancement of the g- 
factor is weak and practically independent of the elec- 
tron density, while the cyclotron mass becomes strongly 
enhanced as the density is decreased. The obtained data 
agree well with the g-factor and effective mass obtained 
by transport measurements, as well as with the Pauli spin 
susceptibility obtained by magnetization measurements 
in parallel magnetic fields, even though the lowest elec- 
tron densities reached in the experiment are somewhat 
higher. Thus, we arrive at the conclusion that, unlike in 
the Stoner scenario, it is indeed the effective mass that 
is responsible for the dramatically enhanced spin suscep- 
tibility at low electron densities. 
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Note added in proof. — After this work had been com- 
pleted, we learned that Punnoose and Finkelstein ^| 
made a renormalization group analysis for multi-valley 
2D systems. Their conclusion that the effective mass 
dramatically increases at the metal-insulator transition 
while the g-factor remains nearly intact is consistent with 
our experimental results. 
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